A Summary of Two Decades of QTL and Candidate Genes That Control Seed Tocopherol Contents in Maize (Zea mays L.).

Tocopherols are secondary metabolites synthesized through the shikimate biosynthetic pathway in the plastids of most plants. It is well known that α-Tocopherol (vitamin E) has many health benefits for humans and animals; therefore, it is highly used in human and animal diets. Tocopherols vary considerably in most crop (and plant) species and within cultivars of the same species depending on environmental and growth conditions; tocopherol content is a polygenic, complex traits, and its inheritance is poorly understood. The objective of this review paper was to summarize all identified quantitative trait loci (QTL) that control seed tocopherols and related contents identified in maize (Zea mays) during the past two decades (2002-2022). Candidate genes identified within these QTL regions are also discussed. The QTL described here, and candidate genes identified within these genomic regions could be used in breeding programs to develop maize cultivars with high, beneficial levels of seed tocopherol contents.

Genetic Mapping for QTL Associated with Seed Nickel and Molybdenum Accumulation in the Soybean 'Forrest' by 'Williams 82' RIL Population.

Understanding the genetic basis of seed Ni and Mo is essential. Since soybean is a major crop in the world and a major source for nutrients, including Ni and Mo, the objective of the current research was to map genetic regions (quantitative trait loci, QTL) linked to Ni and Mo concentrations in soybean seed. A recombinant inbred line (RIL) population was derived from a cross between 'Forrest' and 'Williams 82' (F × W82). A total of 306 lines was used for genotyping using 5405 single nucleotides polymorphism (SNP) markers using Infinium SNP6K BeadChips. A two-year experiment was conducted and included the parents and the RIL population. One experiment was conducted in 2018 in North Carolina (NC), and the second experiment was conducted in Illinois in 2020 (IL). Logarithm of the odds (LOD) of ≥2.5 was set as a threshold to report identified QTL using the composite interval mapping (CIM) method. A wide range of Ni and Mo concentrations among RILs was observed. A total of four QTL (qNi-01, qNi-02, and qNi-03 on Chr 2, 8, and 9, respectively, in 2018, and qNi-01 on Chr 20 in 2020) was identified for seed Ni. All these QTL were significantly (LOD threshold > 2.5) associated with seed Ni, with LOD scores ranging between 2.71-3.44, and with phenotypic variance ranging from 4.48-6.97%. A total of three QTL for Mo (qMo-01, qMo-02, and qMo-03 on Chr 1, 3, 17, respectively) was identified in 2018, and four QTL (qMo-01, qMo-02, qMo-03, and qMo-04, on Chr 5, 11, 14, and 16, respectively) were identified in 2020. Some of the current QTL had high LOD and significantly contributed to the phenotypic variance for the trait. For example, in 2018, Mo QTL qMo-01 on Chr 1 had LOD of 7.8, explaining a phenotypic variance of 41.17%, and qMo-03 on Chr 17 had LOD of 5.33, with phenotypic variance explained of 41.49%. In addition, one Mo QTL (qMo-03 on Chr 14) had LOD of 9.77, explaining 51.57% of phenotypic variance related to the trait, and another Mo QTL (qMo-04 on Chr 16) had LOD of 7.62 and explained 49.95% of phenotypic variance. None of the QTL identified here were identified twice across locations/years. Based on a search of the available literature and of SoyBase, the four QTL for Ni, identified on Chr 2, 8, 9, and 20, and the five QTL associated with Mo, identified on Chr 1, 17, 11, 14, and 16, are novel and not previously reported. This research contributes new insights into the genetic mapping of Ni and Mo, and provides valuable QTL and molecular markers that can potentially assist in selecting Ni and Mo levels in soybean seeds.

Quantitative Trait Loci and Candidate Genes That Control Seed Sugars Contents in the Soybean 'Forrest' by 'Williams 82' Recombinant Inbred Line Population.

Soybean seed sugars are among the most abundant beneficial compounds for human and animal consumption in soybean seeds. Higher seed sugars such as sucrose are desirable as they contribute to taste and flavor in soy-based food. Therefore, the objectives of this study were to use the 'Forrest' by 'Williams 82' (F × W82) recombinant inbred line (RIL) soybean population (n = 309) to identify quantitative trait loci (QTLs) and candidate genes that control seed sugar (sucrose, stachyose, and raffinose) contents in two environments (North Carolina and Illinois) over two years (2018 and 2020). A total of 26 QTLs that control seed sugar contents were identified and mapped on 16 soybean chromosomes (chrs.). Interestingly, five QTL regions were identified in both locations, Illinois and North Carolina, in this study on chrs. 2, 5, 13, 17, and 20. Amongst 57 candidate genes identified in this study, 16 were located within 10 Megabase (MB) of the identified QTLs. Amongst them, a cluster of four genes involved in the sugars' pathway was collocated within 6 MB of two QTLs that were detected in this study on chr. 17. Further functional validation of the identified genes could be beneficial in breeding programs to produce soybean lines with high beneficial sucrose and low raffinose family oligosaccharides.

Soybean gene co-expression network analysis identifies two co-regulated gene modules associated with nodule formation and development.

Gene co-expression network analysis is an efficient systems biology approach for the discovery of novel gene functions and trait-associated gene modules. To identify clusters of functionally related genes involved in soybean nodule formation and development, we performed a weighted gene co-expression network analysis. Two nodule-specific modules (NSM-1 and NSM-2, containing 304 and 203 genes, respectively) were identified. The NSM-1 gene promoters were significantly enriched in cis-binding elements for ERF, MYB, and C2H2-type zinc transcription factors, whereas NSM-2 gene promoters were enriched in cis-binding elements for TCP, bZIP, and bHLH transcription factors, suggesting a role of these regulatory factors in the transcriptional activation of nodule co-expressed genes. The co-expressed gene modules included genes with potential novel roles in nodulation, including those involved in xylem development, transmembrane transport, the ethylene signalling pathway, cytoskeleton organization, cytokinesis and regulation of the cell cycle, regulation of meristem initiation and growth, transcriptional regulation, DNA methylation, and histone modifications. Functional analysis of two co-expressed genes using TILLING mutants provided novel insight into the involvement of unsaturated fatty acid biosynthesis and folate metabolism in nodule formation and development. The identified gene co-expression modules provide valuable resources for further functional genomics studies to dissect the genetic basis of nodule formation and development in soybean.

Proteomic, Transcriptomic, Mutational, and Functional Assays Reveal the Involvement of Both THF and PLP Sites at the GmSHMT08 in Resistance to Soybean Cyst Nematode.

The serine hydroxymethyltransferase (SHMT; E.C. 2.1.2.1) is involved in the interconversion of serine/glycine and tetrahydrofolate (THF)/5,10-methylene THF, playing a key role in one-carbon metabolism, the de novo purine pathway, cellular methylation reactions, redox homeostasis maintenance, and methionine and thymidylate synthesis. GmSHMT08 is the soybean gene underlying soybean cyst nematode (SCN) resistance at the Rhg4 locus. GmSHMT08 protein contains four tetrahydrofolate (THF) cofactor binding sites (L129, L135, F284, N374) and six pyridoxal phosphate (PLP) cofactor binding/catalysis sites (Y59, G106, G107, H134, S190A, H218). In the current study, proteomic analysis of a data set of protein complex immunoprecipitated using GmSHMT08 antibodies under SCN infected soybean roots reveals the presence of enriched pathways that mainly use glycine/serine as a substrate (glyoxylate cycle, redox homeostasis, glycolysis, and heme biosynthesis). Root and leaf transcriptomic analysis of differentially expressed genes under SCN infection supported the proteomic data, pointing directly to the involvement of the interconversion reaction carried out by the serine hydroxymethyltransferase enzyme. Direct site mutagenesis revealed that all mutated THF and PLP sites at the GmSHMT08 resulted in increased SCN resistance. We have shown the involvement of PLP sites in SCN resistance. Specially, the effect of the two Y59 and S190 PLP sites was more drastic than the tested THF sites. This unprecedented finding will help us to identify the biological outcomes of THF and PLP residues at the GmSHMT08 and to understand SCN resistance mechanisms.

QTL and Candidate Genes for Seed Tocopherol Content in 'Forrest' by 'Williams 82' Recombinant Inbred Line (RIL) Population of Soybean.

Soybean seeds are rich in secondary metabolites which are beneficial for human health, including tocopherols. Tocopherols play an important role in human and animal nutrition thanks to their antioxidant activity. In this study, the 'Forrest' by 'Williams 82' (F×W82) recombinant inbred line (RIL) population (n = 306) was used to map quantitative trait loci (QTL) for seed α-tocopherol, ß-tocopherol, δ -tocopherol, γ-tocopherol, and total tocopherol contents in Carbondale, IL over two years. Also, the identification of the candidate genes involved in soybean tocopherols biosynthetic pathway was performed. A total of 32 QTL controlling various seed tocopherol contents have been identified and mapped on Chrs. 1, 2, 5, 6, 7, 8, 9, 10, 12, 13, 16, 17, and 20. One major and novel QTL was identified on Chr. 6 with an R2 of 27.8, 9.9, and 6.9 for δ-tocopherol, α-tocopherol, and total tocopherol content, respectively. Reverse BLAST analysis of the genes that were identified in Arabidopsis allowed the identification of 37 genes involved in soybean tocopherol pathway, among which 11 were located close to the identified QTLs. The tocopherol cyclase gene (TC) Glyma.06G084100 is located close to the QTLs controlling δ-tocopherol (R2 = 27.8), α-tocopherol (R2 = 9.96), and total-tocopherol (R2 = 6.95). The geranylgeranyl diphosphate reductase (GGDR) Glyma.05G026200 gene is located close to a QTL controlling total tocopherol content in soybean (R2 = 4.42). The two methylphytylbenzoquinol methyltransferase (MPBQ-MT) candidate genes Glyma.02G002000 and Glyma.02G143700 are located close to a QTL controlling δ-tocopherol content (R2 = 3.57). The two γ-tocopherol methyltransferase (γ-TMT) genes, Glyma.12G014200 and Glyma.12G014300, are located close to QTLs controlling (γ+ß) tocopherol content (R2 = 8.86) and total tocopherol (R2 = 5.94). The identified tocopherol seed QTLs and candidate genes will be beneficial in breeding programs to develop soybean cultivars with high tocopherol contents.

The Soybean High Density 'Forrest' by 'Williams 82' SNP-Based Genetic Linkage Map Identifies QTL and Candidate Genes for Seed Isoflavone Content.

Isoflavones are secondary metabolites that are abundant in soybean and other legume seeds providing health and nutrition benefits for both humans and animals. The objectives of this study were to construct a single nucleotide polymorphism (SNP)-based genetic linkage map using the 'Forrest' by 'Williams 82' (F×W82) recombinant inbred line (RIL) population (n = 306); map quantitative trait loci (QTL) for seed daidzein, genistein, glycitein, and total isoflavone contents in two environments over two years (NC-2018 and IL-2020); identify candidate genes for seed isoflavone. The FXW82 SNP-based map was composed of 2075 SNPs and covered 4029.9 cM. A total of 27 QTL that control various seed isoflavone traits have been identified and mapped on chromosomes (Chrs.) 2, 4, 5, 6, 10, 12, 15, 19, and 20 in both NC-2018 (13 QTL) and IL-2020 (14 QTL). The six QTL regions on Chrs. 2, 4, 5, 12, 15, and 19 are novel regions while the other 21 QTL have been identified by other studies using different biparental mapping populations or genome-wide association studies (GWAS). A total of 130 candidate genes involved in isoflavone biosynthetic pathways have been identified on all 20 Chrs. And among them 16 have been identified and located within or close to the QTL identified in this study. Moreover, transcripts from four genes (Glyma.10G058200, Glyma.06G143000, Glyma.06G137100, and Glyma.06G137300) were highly abundant in Forrest and Williams 82 seeds. The identified QTL and four candidate genes will be useful in breeding programs to develop soybean cultivars with high beneficial isoflavone contents.

Genome-wide identification and analysis of soybean acyl-ACP thioesterase gene family reveals the role of GmFAT to improve fatty acid composition in soybean seed.

KEY MESSAGE: Soybean acyl-ACP thioesterase gene family have been characterized; GmFATA1A mutants were discovered to confer high oleic acid, while GmFATB mutants presented low palmitic and high oleic acid seed content. Soybean oil stability and quality are primarily determined by the relative proportions of saturated versus unsaturated fatty acids. Commodity soybean typically contains 11% palmitic acid, as the primary saturated fatty acids. Reducing palmitic acid content is the principal approach to minimize the levels of saturated fatty acids in soybean. Though high palmitic acid enhances oxidative stability of soybean oil, it is negatively correlated with oil and oleic acid content and can cause coronary heart diseases for humans. For plants, acyl-acyl carrier protein (ACP) thioesterases (TEs) are a group of enzymes to hydrolyze acyl group and release free fatty acid from plastid. Among them, GmFATB1A has become the main target to genetically reduce the palmitic acid content in soybean. However, the role of members in soybean acyl-ACP thioesterase gene family is largely unknown. In this study, we characterized two classes of TEs, GmFATA, and GmFATB in soybean. We also denominated two GmFATA members and discovered six additional members that belong to GmFATB gene family through phylogenetic, syntenic, and in silico analysis. Using TILLING-by-Sequencing+, we identified an allelic series of mutations in five soybean acyl-ACP thioesterase genes, including GmFATA1A, GmFATB1A, GmFATB1B, GmFATB2A, and GmFATB2B. Additionally, we discovered mutations at GmFATA1A to confer high oleic acid (up to 34.5%) content, while mutations at GmFATB presented low palmitic acid (as low as 5.6%) and high oleic acid (up to 36.5%) phenotypes. The obtained soybean mutants with altered fatty acid content can be used in soybean breeding program for improving soybean oil composition traits.

TILLING-by-Sequencing+ Reveals the Role of Novel Fatty Acid Desaturases (GmFAD2-2s) in Increasing Soybean Seed Oleic Acid Content.

Soybean is the second largest source of oil worldwide. Developing soybean varieties with high levels of oleic acid is a primary goal of the soybean breeders and industry. Edible oils containing high level of oleic acid and low level of linoleic acid are considered with higher oxidative stability and can be used as a natural antioxidant in food stability. All developed high oleic acid soybeans carry two alleles; GmFAD2-1A and GmFAD2-1B. However, when planted in cold soil, a possible reduction in seed germination was reported when high seed oleic acid derived from GmFAD2-1 alleles were used. Besides the soybean fatty acid desaturase (GmFAD2-1) subfamily, the GmFAD2-2 subfamily is composed of five members, including GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E. Segmental duplication of GmFAD2-1A/GmFAD2-1B, GmFAD2-2A/GmFAD2-2C, GmFAD2-2A/GmFAD2-2D, and GmFAD2-2D/GmFAD2-2C have occurred about 10.65, 27.04, 100.81, and 106.55 Mya, respectively. Using TILLING-by-Sequencing+ technology, we successfully identified 12, 8, 10, 9, and 19 EMS mutants at the GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E genes, respectively. Functional analyses of newly identified mutants revealed unprecedented role of the five GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E members in controlling the seed oleic acid content. Most importantly, unlike GmFAD2-1 members, subcellular localization revealed that members of the GmFAD2-2 subfamily showed a cytoplasmic localization, which may suggest the presence of an alternative fatty acid desaturase pathway in soybean for converting oleic acid content without substantially altering the traditional plastidial/ER fatty acid production.

TILLING-by-Sequencing+ to Decipher Oil Biosynthesis Pathway in Soybeans: A New and Effective Platform for High-Throughput Gene Functional Analysis.

Reverse genetic approaches have been widely applied to study gene function in crop species; however, these techniques, including gel-based TILLING, present low efficiency to characterize genes in soybeans due to genome complexity, gene duplication, and the presence of multiple gene family members that share high homology in their DNA sequence. Chemical mutagenesis emerges as a genetically modified-free strategy to produce large-scale soybean mutants for economically important traits improvement. The current study uses an optimized high-throughput TILLING by target capture sequencing technology, or TILLING-by-Sequencing+ (TbyS+), coupled with universal bioinformatic tools to identify population-wide mutations in soybeans. Four ethyl methanesulfonate mutagenized populations (4032 mutant families) have been screened for the presence of induced mutations in targeted genes. The mutation types and effects have been characterized for a total of 138 soybean genes involved in soybean seed composition, disease resistance, and many other quality traits. To test the efficiency of TbyS+ in complex genomes, we used soybeans as a model with a focus on three desaturase gene families, GmSACPD, GmFAD2, and GmFAD3, that are involved in the soybean fatty acid biosynthesis pathway. We successfully isolated mutants from all the six gene family members. Unsurprisingly, most of the characterized mutants showed significant changes either in their stearic, oleic, or linolenic acids. By using TbyS+, we discovered novel sources of soybean oil traits, including high saturated and monosaturated fatty acids in addition to low polyunsaturated fatty acid contents. This technology provides an unprecedented platform for highly effective screening of polyploid mutant populations and functional gene analysis. The obtained soybean mutants from this study can be used in subsequent soybean breeding programs for improved oil composition traits.

Soybean TILLING-by-Sequencing+ reveals the role of novel GmSACPD members in unsaturated fatty acid biosynthesis while maintaining healthy nodules.

Developing soybean lines with high levels of stearic acid is a primary goal of the soybean industry. Most high-stearic-acid soybeans carry different GmSACPD-C mutated alleles. However, due to the dual role of GmSACPD-C in seeds and nodule development, all derived deleterious GmSACPD-C mutant alleles are of extremely poor agronomic value because of defective nodulation. The soybean stearoyl-acyl carrier protein desaturase (GmSACPD) gene family is composed of five members. Comparative genomics analysis indicated that SACPD genes were duplicated and derived from a common ancestor that is still present in chlorophytic algae. Synteny analysis showed the presence of segment duplications between GmSACPD-A/GmSACPD-B, and GmSACPD-C/GmSACPD-D. GmSACPD-E was not contained in any duplicated segment and may be the result of tandem duplication. We developed a TILLING by Target Capture Sequencing (Tilling-by-Sequencing+) technology, a versatile extension of the conventional TILLING by sequencing, and successfully identified 12, 14, and 18 ethyl methanesulfonate mutants at the GmSACPD-A, GmSACPD-B, and GmSACPD-D genes, respectively. Functional analysis of all identified mutants revealed an unprecedented role of GmSACPD-A, GmSACPD-B, and GmSACPD-D in unsaturated fatty acid biosynthesis without affecting nodule development and structure. This discovery will positively impact the development of high-stearic-acid lines to enhance soybean nutritional value without potential developmental tradeoffs.

Mutations at the Serine Hydroxymethyltransferase Impact its Interaction with a Soluble NSF Attachment Protein and a Pathogenesis-Related Protein in Soybean.

Resistance to soybean cyst nematodes (SCN) in "Peking-type" resistance is bigenic, requiring Rhg4-a and rhg1-a. Rhg4-a encodes a serine hydroxymethyltransferase (GmSHMT08) and rhg1-a encodes a soluble NSF attachment protein (GmSNAP18). Recently, it has been shown that a pathogenesis-related protein, GmPR08-Bet VI, potentiates the interaction between GmSHMT08 and GmSNAP18. Mutational analysis using spontaneously occurring and ethyl methanesulfonate (EMS)-induced mutations was carried out to increase our knowledge of the interacting GmSHMT08/GmSNAP18/GmPR08-Bet VI multi-protein complex. Mutations affecting the GmSHMT08 protein structure (dimerization and tetramerization) and interaction sites with GmSNAP18 and GmPR08-Bet VI proteins were found to impact the multi-protein complex. Interestingly, mutations affecting the PLP/THF substrate binding and catalysis did not affect the multi-protein complex, although they resulted in increased susceptibility to SCN. Most importantly, GmSHMT08 and GmSNAP18 from PI88788 were shown to interact within the cell, being potentiated in the presence of GmPR08-Bet VI. In addition, we have shown the presence of incompatibility between the GmSNAP18 (rhg1-b) of PI88788 and GmSHMT08 (Rhg4-a) from Peking. Components of the reactive oxygen species (ROS) pathway were shown to be induced in the SCN incompatible reaction and were mapped to QTLs for resistance to SCN using different mapping populations.